The cutoff-dependence of the Casimir force within an inhomogeneous medium 
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We consider the ground state energy of the electromagnetic field in a piston geometry. In the 
idealised case, where the piston and the walls of the chamber are taken as perfect mirrors, the 
Casimir pressure on the piston is finite and independent of the small scale physics of the media that 
compose the mirrors; the Casimir-energy of the system can be regularised and is cutoff-independent. 
Yet we find that, when the body of the piston is filled with an inhomogeneous dielectric medium, 
the Casimir energy is cutoff-dependent, and the value of the pressure is thus inextricably dependent 
on the detailed behaviour of the mirror and the medium at large wave-vectors. This result is 
inconsistent with recent proposals for regularising Casimir forces in inhomogeneous media. 

PACS numbers: 42.50.-p,42.50.Lc 
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I. INTRODUCTION 

It is now a well-known and experimentally supported 
fact that two parallel uncharged mirrors at zero temper- 
ature will exert an attractive force upon each other, as 
Casimir predicted [T] — a force arising from the ground 
state properties of the electromagnetic field. The the- 
ory has been made quantitative, applying to media de- 
scribed by e(w) and /i(a;) satisfying the Kramers-Kronig 
relations |2H5]: and Casimir forces have been calculated 
for a variety of systems and geometries (e.g. see [SI [7]). 
The purpose of the following discussion is to demonstrate 
that there remains a problem within the theory of the 
Casimir effect for the case of inhomogeneous media, in 
spite of recent efforts to solve or circumvent it [BJ ISj • In 
a separate paper |10j we have shown that this problem is 
inherent within the Lifshitz theory of the Casimir effect, 
where dispersion and dissipation are properly included; 
the Casimir-stress in such cases is infinite and resists reg- 
ularisation. 

But this problem is not peculiar to Lifshitz theory. 
Here we show that, even when wc make the most basic 
attempt to mimic dielectric media with non-dispersive 
boundary conditions — considering only a simple energy 
summation of the field modes — the Casimir force is infi- 
nite in the limit where the regularization (cut-off) tends 
to infinite frequency. It seems apparent that we cannot 
obtain an expression for the zero point force that is inde- 
pendent of the choice of cut off in the energy summation. 



II. THE CASIMIR PISTON 

We first briefly illustrate the principles of our main cal- 
culation in an idealized system where we know there ex- 
ists a finite expression for the Casimir force, independent 
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of the microscopic physics of the bodies involved. The 
chosen system is a cavity of length L, and cross-sectional 
area A = LyLz, divided by a mirror at a distance a from 
the left-most cavity wall (figure [T]). Even at T = OK the 
moveable mirror situated at a; = a is subject to a force 
due to the dependence of the field energy on a. 

We can formally write the ground state energy of the 
system as the sum of all possible contributions of tujj/2, 
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where A G {1, 2} indicates the polarization, and the range 
of summation for each of the indices runs over the allowed 
modes. The eigenfrequencies of the cavities to the left {L) 
and right (i?) of a are independent of the polarization 
when m > 0, 
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and when m = 0, the A = 1 polarization is not an allowed 
mode of the cavity. Given the behaviour of ([2]), ([T]) is 
not a meaningful expression: the summand becomes ever 
larger as m, p & g are increased, and the sum diverges. 
This divergence is due to the artificial assumption that 
there exist mirrors that act at all frequencies. To fix this, 
a factor is inserted into the eigenfrequencies ^ to make 
the sum converge. 
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where f is a free parameter, and a factor of two comes 
from the symmetry of the system with respect to polar- 
ization. This modified expression for the energy can evi- 
dently no longer be considered as the total energy of the 
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system of field plus mirrors. Taken literally as such, ([3| 
would imply that the eigenfrequencies of the field eventu- 
ally all tend to zero, an assumption for which there is no 
obvious motivation. Instead, ([s]) should be interpreted 
as the part of the total energy ([T]) associated with the 
configuration of mirrors in the cavity (the energy avail- 
able to do work on the mirrors) . The exponential factors 
then amount to a model for the dispersive behaviour of 
the mirrors: i.e. at sufficiently high frequencies/ wave- 
vectors the mirrors become transparent, and after this 
the total energy does not depend on their configuration. 
Note that the factor of two in (|3| takes into account the 
sum over polarization. This is not correct for the m = 
mode. However, the energy of this mode is independent 
of a. 



Q we apply the identity. 
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yielding, 
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(6) 

The summation within ([6| is evidently a geometric se- 
ries, which can be evaluated, J2m=o^~"™ ~ 1/(1 
giving, 




FIG. 1: Schematic of the Casimir piston. Two fixed mirrors 
are positioned at a; = and x = L, enclosed by reflecting 
walls a.t y — ±Ly/2 and z — ±Lz/2 (dashed lines). Within 
the chamber is vacuum, and a moveable mirror at a; = a. 
While the ground state energy of the electromagnetic field in 
this system depends on the small scale physics of the mirrors, 
the part of it dependending on a does not. 



To explicitly evaluate ([3| we take the limit 
Ly /a, Lz/a ^ oo, where the summation over m and p can 
be converted into an integration over ky and kz, A.ki — 
t: / Li — > dki (i is either y ov z). After an integration over 
the angle, of fey = kyij -f k^z = k\\ [cos(f?)y -I- sin(6')2;], 
the energy per unit area of the piston configuration is 
found to be 
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The introduction of the exponential factor into Q rep- 
resents an extremely artificial model for the behaviour of 
the mirrors at high frequencies. We therefore separate 
the energy into those parts that depend on ^, and those 
that do not. In anticipation of taking the limit ^ — >■ 0, 
the quantity to the right of the square brackets in ([t]) is 
expanded as far as 
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Inserting Q into Q gives finally, 
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As expected, the energy becomes increasingly large as 
^ — 0. However, the part of the energy depending on 
the position of the mirror is independent of ^. One may 
interpret this to mean that so long as ^/a is negligibly 
small (i.e. neglecting the positive powers of ^ in ^ is 
legitimate), the part of the energy that depends on a is 
independent of how the mirror becomes transparent at 
high frequencies/ wave-vectors. Taking the derivative of 
([9]) with respect to a then yields the usual finite and ^ 
independent pressure. 
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so that the Casimir force on the mirror pulls it towards 
the closer of the two end walls of the chamber. 



III. THE INHOMOGENEOUS CASIMIR 
PISTON 



where a — ^ i — a indicates a repetition of the previous 
expression with 'a' replaced everywhere by 'L — a'. To 



The above procedure for making ([T]) convergent yields 
terms dependent on a that are all finite or zero in the 
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limit ^ — >■ 0. Yet it is not obvious that this fortuitous 
situation occurs for fundamental reasons. Suppose we 
take the same cavity with an inhomogeneous medium 
within the chamber: are the divergent terms in E still 
independent of a? 

Within the chamber we assume a permeability and per- 
mittivity given by (e.g. see figure [2|, 

/i(a;) = ^0 

e{x) = eo [1 + 5e{x)] . (10) 

It is imagined that the medium with permittivity and 
permeability given by ( 10 ) is a rigid body. Is it possible to 



express the force holding the mirror fixed at a (opposing 
the Casimir force) in a form that depends only on e and 
/i? 




the region x £ [0, a] are given by, 
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and, 
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where we assume the limit of Ly/a,Lz/a — >■ oo as in 
section [llj After substituting, a ^ L ~ a, the modes in 
the region x a [a, L] are also given by ( 12 13 ). All modes 
are normalised over the volume, V ^ of each region of the 
chamber, \E^^J^(fx = 1. 

We consider the case where 5e{x) ^ 1, and write the 
frequency of each eigenmode as uim,\ = ^. 
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'Jrn,\ = ^m.\ ^ -^m,A' whcrc the quantities 
containing a superscript '(1)' are supposed linear in 5e. 
To first order in (5e, ( 11 ) is 



the field as 
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FIG. 2: As in figure[T] we have a perfectly reflecting rectangu- 
lar chamber of length L and cross sectional area, A. Within 
the chamber is a mirror at a; = a, surrounded by an inhomo- 
geneous dielectric with e and /i given by |10 1. We look for the 



dependence of the energy of this system on a. 

As in the previous section, and in common with 
Casimir's original calculation [T], we make the artificial 
assumption that e and are independent of frequency. 
This needs some justification. These effects are of course 
fundamental to the interaction of light and matter |llj . 
To properly account for them requires Lifshitz theory [21- 
[5, and in separate work we have shown that finite an- 
swers cannot be obtained from this theory in the above 
case [ini. The purpose of this calculation is to exam- 
ine whether this problem is peculiar to Lifshitz theory, 
or still present in a much more naive theory where the 
propagation speed of light is simply varied from point to 
point, and not as a function of frequency. 

The situation is slightly complicated in comparison to 
section [n] due to the fact that the two polarizations do 
not behave degenerately in the medium defined by (10 1. 
Within the cavity, the electromagnetic field obeys. 



[1 + Se{x)]Ern,x ^ 0, (11) 



Multiplying ( 14 ) on the left by -E^*^"*^ and integrating over 



V (which couM either be the left or right region of fig- 
ure [2|) after an integration by parts we find, 
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which is the standard expression for the first order per- 
turbation of the eigenfrequencies of an optical cavity (see 
e.g. [I2t 1131). For the particular case, 
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we find after inserting ([12j-[13j) into ( |15[ ), that in the left 
most portion of the piston with A = 1, 
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where m labels the spatial dependence of the mode, and 
A € {1,2} the polarization. When Se = 0, the modes in 



after substituting a ^ L — a into (17 18), one has the 



expressions for the right hand side of the cavity. The 
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above shifts in the eigenfrequencies are independent of 
the angle, 9 of fcy. We can therefore perform the inte- 
gration over 9 in the mode summation. The change in 
energy per unit area due to the presence of the medium 
within the piston is computed, and takes the form, 



unperturbed case. The integration over fcy is performed 
using the same techniques as in section [llj The results are 
given by ( A2 - A3 1 in appendix ^ In addition to terms 
proportional to geometric series"c.f. ([6])), the resulting 
sum over m contains terms of the form, 



A Ait 



k\\dk\\ 



.Ml),-,.--/. 



R{1) 



(19) 



where we use the same regularizing function as in the 



3>(e-«-/^l,^;)= ^ 



m -\- V 



(20) 



The quantity $ is known as the Lerch function [T31 US] . 
Applying the notation of (20), we can perform the sum- 
mation, with the result, 
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where we have converted the approximate expression. 



( 19 1 into an exact relation through taking the limit 
a ^ da (as is done e.g. in [E]). Equation (21| rep- 



resents the rate of change of E as the amplitude of ( 16 ) 



is increased from zero to da. Again we consider the limit 
f — 0, and make use of the following series expansion of 
the Lerch function (18j . 

$(e-«''/", 1, v)e-^''"/'' + log(C7r/a) = -7 - ^p{v) 

m+l 



' m + l (m + l)! ^ ^ 

m— 



where 7 is Euler's constant, 'ip{a/2L) the digamma func- 
tion (logarithmic derivative of the gamma function) , and 



the Bm+i{a/2L) are Bernoulli polynomials [T3]. From 
(22 1 we find, after neglecting positive powers of ^, 
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Applying (23 1 and (|8|) to (21 1, gives the following lengthy 



expression for the rate of change of energy with respect 
to a, 
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Evidently the divergent terms — except those propor- 
tional to and — are dependent on the position 
of the piston within the chamber. These divergences are 
not cancelled by including the contributions to the energy 
from the other side of the piston. This means that in the 
limit ^ — >■ 0, the force on the piston is discontinuous as 
a function of a, being finite when a = 0, and infinite as 



a is moved away from zero. It therefore seems that a 
^ independent meaning cannot be given to the Casimir 
force, when in an inhomogeneous medium. The value of 
the force is cutoff-dependent. 
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IV. RELATION TO PREVIOUS PROPOSALS 

It was recently proposed [9] that one should calculate 
the Casimir force in an inhomogeneous medium through 
forming a Laurent expansion of the energy in powers of 
the regularizing parameter (here ^). The regularized en- 
ergy is then defined to be "t/ie term, cq in [the] Laurent 
expansion that is independent of [^] and corresponds to 
discarding the principal part of the Laurent series before 
taking the limit as [^] tends to zerd'^ . This procedure 
works well for homogeneous systems. However, apply- 
ing this procedure to ( [24| (which is proportional to the 
change in energy due to the inhomogeneous medium, 
when a is small), proves problematic. For one thing, 
we would still be left with the logarithmic divergence. 
More seriously, however, the negative powers of ^ which 
depend on a are simply swept away. Therefore the reg- 
ularization procedure modifies the value of the force on 
the mirrors. There is no obvious physical justification for 
this and we emphasize, as we did in section [llj that the 
regularization in the case of a homogeneous medium does 
not modify the value of the force (the diverging terms are 
independent of the positions of the mirrors) . 

V. CONCLUSIONS 

It is has been demonstrated in previous work |10] 
that there is a problem with the Lifshitz theory of 



Casimir forces: the Casimir-stress in inhomogeneous me- 
dia, where the optical properties of the medium vary con- 
tinuously in space, is infinite and resists regularisation. 
But the problem of computing Casimir forces in inho- 
mogeneous media appears to be systemic. Here wc have 
shown that, even with a simple energy mode summa- 
tion, the ground-state energy of a system proves simi- 
larly resistive to regularisation. We have shown this for 
the case of a Casimir piston, when an inhomogeneous 
medium is present in the cavity. Our calculation demon- 
strates a cutoff dependence in the Casimir force which 
suggests, surprisingly, that the Casimir forces in a system 
depend in detail on its microphysical properties. If this 
is the case, it seems unlikely that a generally finite and 
physically meaningful result could be obtained through 
a simple modification to the existing regularisation pro- 
cedure [9l [17] ; some additional physics must be taken 
into account. 
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Appendix A: Integration of perturbed eigenfrequencies 

Here we give the results of the following integrals. 
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18), with the s — 'i?' expressions are obtained after the substitution a ^ L 



The eigenfrequency shifts for polarization, A = 1 depends on only through uj'^^i'. i.e. in the same way as the 
unperturbed eigenfrequencies. Therefore (Al) can be performed for this polarization in the same way as section [lH 
with the result. 
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Polarization, A = 2, has a slightly more complicated dependence on fcy. Yet there is nothing fundamentally different 
about performing the integrals, and one obtains. 
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After substituting a 



L — a into (A2) and (A3), one obtains ^ and L^ 2- 
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